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Abstract 

We review inflationary cosmology in modified gravity such as R 2 gravity with its extensions in 
order to generalize the Starobinsky inflation model. In particular, we explore inflation realized 
by three kinds of effects: modification of gravity, the quantum anomaly, and the R 2 term in loop 
quantum cosmology. It is explicitly demonstrated that in these inflationary models, the spectral 
index of scalar modes of the density perturbations and the tensor-to-scalar ratio can be consistent 
with the Planck results. Bounce cosmology in F(R) gravity is also explained. 
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I. INTRODUCTION 


Inflation in the early universe has recently been studied much more extensively because 


of the BICEP2 experiment [l| in terms of the primordial gravitational waves, in addition 


to the Wilkinson Microwave anisotropy probe (WMAP) 2-6] and the Planck satellite 0,0 
on the unisotropy of the cosmic microwave background (CMB) radiation. For a standard 
inflationary scenario like chaotic inflation jfi], the existence of the inflaton field is assumed, 
whose potential contributes to inflation. 

On the other hand, the accelerated expansion of the universe including inflation and the 
late-time acceleration, i.e., dark energy problem, can be realized in modified gravity theories 
such as F(R) gravity (for reviews on inflation, see, e.g., {lO, 11], whereas for dark energy and 


modified gravity, see, for example, Refs 
inflation such as the Starobinsky inflation 


21 


2Qf 


For instance, the trace-anomaly driven 


22] is well known. 


Indeed, the WMAP and Planck data [2-8] support a kind of the trace-anomaly driven 
inflation with the R 2 term. Such a theory can be regarded as modified gravity because the 
R 2 term or its higher derivative term of the trace-anomaly term QR, which leads to the long 
enough inflation and graceful exit from it [21|, is the effective action of gravity, where □ is 
the covariant d’Alembertian for a scalar quantity 1 . 


In this paper, we review the main results in Refs. 


26 


The main purpose of this 


paper is to explain the recent developments on inflationary models to realize the Planck 
results in the so-called R 2 gravity (namely, the action consist of the Einstein-Hilbert term 
plus R 2 term) with further extensions, which can be regarded as a kind of F(R) gravity. 
Particularly, we consider inflation (i) derived by modification terms of gravity 26], (ii) 


through the quantum anomaly 


27L and (iii) in R 2 gravity in the framework of the so-called 


loop quantum cosmo 


see, for example, 


34 


ogy (LQC) [29M33|] to include quantum effects [28] (for reviews on LQC, 


40|). In addition, we state the recent progress of the bounce cosmology 


in F(R ) gravity by presenting the important consequences in Refs. 4ll, [42|. We use units 
of k-Q = c — h — 1 and express the gravitational constant 87 tGn by k, 2 = 8tt/Mp\ 2 with the 
Planck mass of Mpi = = 1.2 x 10 19 GeV. 

The organization of the paper is the following. In Sec. II, we study inflation induced by 


23-25j, the features of inflation in non-local gravity including such a non-local term as \3R has 


In Refs. 


been analyzed in detail. 
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modification of gravity. In Sec. Ill, we explore the trace-anomaly driven inflation in modified 
gravity. In Sec. IV, we investigate R 2 gravity in the context of LQC. Furthermore, in Sec. V, 
we reconstruct F(R) gravity to realize the cosmological bounce in LQC. Finally, conclusions 
are described in Sec. VI. 


II. INFLATION INDUCED BY MODIFICATION OF GRAVITY 


In this section, we review inflation in modified gravity, particularly F(R) gravity, based 


on Ref. 


26| . The deviation of F(R) gravity from general relativity may be interpreted as a 


kind of quantum corrections in the early universe, or such a modification of gravity could 
be motivated by the so-called ultraviolet (UV) completion of quantum gravity. In fact, the 


Starobinsky inflation 21] can be regarded as inflation induced by the modification term of 
R 2 from general relativity. We here attempt to examine inflation by the other forms of 
modification of gravity. 


A. Conformal transformation 


We first explain the conformal transformation from F(R) 


gravity in the Jordan frame 


to the corresponding scalar held theory in the Einstein frame 19, |43]. The action of F(R) 


gravity is represented as S — J d A x^/—g [F(R)/ (2k 2 )], where g is the determinant of the 
metric tensor g /lu . We use a conformal transformation g^ u = Vt~ 2 g^ v with fl 2 = F r , where 
the hat denotes quantities in the Einstein frame, and the subscription of Fr denotes the 
derivative with respect to R as F r (R) = dF(R)/dR. Here, we introduce a scalar held 
ip = —yQ/2 (Q/k) In Fr. Through the conformal transformation, the action in the Einstein 
frame reads 


44 


45] 


s E = J f 4L _ Ig^a^a, 


,if - V{ip) 


V{<p) = 


FrR - F 


2k 2 (F r ) 


2 > 


(Hi) 

(II.2) 


with F r = exp ^2/3 Kip^j. This is the action for a canonical scalar held ip with its 
potential V{p). For the Starobinsky inhation model 2JL] with F(R) — R + ctgK 2 /? 2 , where 
as is a constant, we have V(<p) = [1/ (8asK 2 )] (1 — exp ( — \/2/2>Kp 
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B. Slow-roll inflation 


We describe the procedures to deal with the so-called slow-roll inflation. We consider 
the action in Eq. (III.ID and regard ip as the inflaton field. We suppose the flat Friedmann- 
Lemaitre-Robertson-Walker (FLRW) metric ds 2 = — dt 2 + a 2 (t)dx 2 . Here, a(t) is the scale 
factor. The Hubble parameter is defined as H = a/a, where the dot denotes the time 
derivative. The gravitational field equations in this background read 3 H 2 / n 2 = ip 2 /2 + V(ip) 
and — (^2H + 3H 2 ^J /k 2 = ip 2 / 2 — V{ip). Moreover, the equation of motion (EoM) for ip 
becomes ip + 3 Hip + dV(ip)/dip = 0. 

For the slow-roll regime, we impose the slow-roll approximations of ip 2 /2 <C V(ip) 
on the Friedmann equation and \ip\ <C \3Hip\ on the EoM for p , so that we can find 
3 H 2 /k 2 « V(ip) ~ constant and 3 Hip + dV (ip) /dip ~ 0. Furthermore, we define the slow-roll 
parameters e = —H/H 2 = [1/ (2k 2 )] [(dV (ip) / dip) /V(ip)} 2 (-C 1) and r/ = —H/ ( 2 HH^j = 
[1/k 2 ] [(d 2 V(ip)/dip 2 ) /V{ip)\ (<S 1). During inflation, these parameters should be much 
smaller than unity. In addition, the number of e-folds is N e = ln(af/a;) = [*/ Hdt ps 
k 2 \V/ (dV(ip) / dip)\ dip. Here, a ; (ip/) and Of (ipt) are the values of the scale factor a (the 
scalar field ip) at the initial time f; and end of time £f of inflation, respectively. Moreover, 
in deriving the second approximate equality, we have used the second gravitational field 
equation with the slow-roll approximation. The amplitude of the power spectrum for the 
curvature perturbations is expressed as = n 2 H 2 / (87r 2 e) ~ k 4 R/ (247r 2 e), where in the 


second approximate equality follows from the Friedmann equation operatec 


proximation, and n s and r are written as n s — 1 = —6e + 2^ and r = 16e 


46 


le slow-roll ap- 


47 ], The detailed 


explanations on the reconstruction of potential of inflationary models has been executed in 
Ref. 48]. 


C. Reconstruction of F(R) gravity 

There are two possible ways to reconstruct F(R) gravity models to realize inflation. One 
is to start from the action in the Einstein frame. The other is to reconstruct the action in 
the Jordan frame. In this subsection, we consider the former way. We also explain the latter 
way in the next subsection. 

The main purpose of our investigations is that we study the generalization of the Starobin- 
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sky inflation model. To execute it, we take an appropriate form of V(<p) in the Jordan 
frame, which is an extended form from that in the Starobinsky inflation model. By taking 
the derivative of Eq. (III. 21) with respect to R, we find 


R.Fji = — \/ 6 k — ( 1 /( 92 ) exp 
dp 



(IL3) 


Through this equation, i.e., Eq. (111.21) . we reconstruct the form of F(R) in the Jordan frame 
from the potential V ( 92 ) in the Einstein frame. 


1. Extension of the Starobinsky inflation model 


As the simplest model, we explore the following potential 


V ( 92 ) = ci + c 2 exp ( J^Kip j + c 3 exp f 2 y ) , 


(IL4) 


where C\ 0), c 2 , and c 3 are constants. In this case, from Eq. (III. 3D we have 2ciF^ + c 2 ib f — 
RFr = 0. By solving this equation with F R 0, we eventually find that the corresponding 
form of F(R) can be expressed as 

r2 

F(R) = R + — + ci - c 3 . (II.5) 

4ci 

Here, we have set —c 2 / (2ci) = 1 in order to reproduce the Einstein-Hilbert term in F(R) 
and used Eq. (III. 2D in determining the integration constant. If C\ = c 3 (which leads to 
c 3 = —c 2 /2), this model is equivalent to the Starobinsky inflation model with V(p) = 
Ci ^1 — exp ^/2/3 k</?^ j , whereas for ci > c 3 , this corresponds to an extended model 
of the Starobinsky inflation with a cosmological constant. As the possibile origins of Such 
a cosmological constant emerging at the large curvature regime could originate from the 
quantum effects, or a modification term of gravity removing the cosmological constant at 


49H52], 


the small curvature regime 

In the following, we set c 3 = 0 for simplicity and introduce a positive 71 (> 0) to express 
Ci as Ci = 71 / (4k 2 ). Here, 71 has the mass dimension 2 and the dimensionless quantity 
71 /Mpi < 1, so that in the higher-curvature regime, the correction to the Einstein gravity 
can appear. We explore the inflationary dynamics in this extended model with the potential 
V{<p) = 71 / (4k 2 ) — [ 71 / (2k 2 )] exp (v^73k^) . We consider the case that the inflaton slowly 
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rolls from the initial value with its large negative amplitude down to the minimum of the 
potential as V(<p = 0) = —71/ (4k 2 ) (< 0). In this case, from the gravitational held equations 
we find that the exponential inflation can be realized as a(t) = a- x exp with the Hubble 

parameter H m f = (1/2) \f^i /3 during infl ation and a; a constant. Furthermore, the solution 
of tp reads ip = — \j3/2 k In (1/3) y/ 271/3 (t; — t ) . Around the beginning of inflation t ~ ti, 

1 — exp (— 

,_ x -1 

1 — exp ( — 


\<p\ 1 and the slow-roll parameters are e = (4/3) 

M = (4/3) 


-2 


<C 1 and 

<C 1. These slow-roll parameters become of order of 
unity when p> approaches ipt ~ —0.17^3/2//?. For |<^i| |y>f|, the number of e-folds is given 

by N e ~ (1/2) a/ 7 'i/ 6 fl. In addition, we find tf — t; — 3^/3/ ( 27 ^ exp (0.17). For fV e = 60, 
we have (p\ ~ 1.07Mpi. The slow-roll parameters are also represented as e ~ 3/ (4A1|) and 
\r]\ ps 1 /A r . As a consequence, we obtain 


A 


k 2 7iA ( 


•ft 


72ft 2 


2 12 

n s — 1 ~ — — , r —— . 




A 2 


(IL6) 


Here, we remark that ~ k 2 7 i A e 2 / (72ft 2 ) -C K 2 M P iA e 2 / (72ft 2 ) because 71 M P \. 


The observations obtained from the Planck satellite suggest n s = 0.9603±0.0073 ( 68 % CL) 


and r < 0.11 (95% CL) [7]. In this model, for n s < 1 and r < 0.11, we see that n s > 
1 — yd).11/3 = 0.809. Accordingly, for N e = 60, we acquire n s = 0.967 and r = 3.00 x 10 3 . 
Thus, in this model, the spectral index n s of the curvature perturbations and the tensor- 
to-scalar ratio r consistent with the Planck result can be realized. Various descriptions of 
inflationary models in terms of scalar field models [53] and perfect fluid as well as F(R) 
gravity [54} have been examined. Moreover, the effects of quantum corrections on inflation 


have been explored in Refs. 


55M57}. We note that the BICEP2 experiment has recently 


detected the R-mode polarization of the cosmic microwave background (CMB) radiation 


with the tensor to scalar ratio r = OAO/qq^ ( 68 % CL) [1]. There have been proposed several 
discussions on the method to obtain this result regarding the subtraction of the foreground 
data, e.g ., Refs. 58H62], A study to support the BICEP2 results has also been reported in 


Ref. 63]. Very recently, the collaboration between BICEP2/Keck and Planck has released 
the result of r < 0.12 for the wave number k = 0.05 Mpc -1 of tensor mode of the density 


perturbations 


64]- 


6 






































2. Power-law corrections to general relativity 


Next, we examine the following potential. 


™Hs) 




1 - exp 1 \I -Kip 


9/(9-l) 


q~ 2 

g- 1 


+i? c exp | \/ -Kip | | exp 


3^ 


exp 


A p exp 2 \ / -Kip 


Kip 


q -1 


(IL7) 


For this potential in the Einstein frame, a model in which a generic power-law correction 
term is added to the Einstein-Hilbert term is reconstructed as 


FIR) ^R+PlR+Wy + Ap- 


(II.8) 


Here, /3 (> 0) is a dimensionful positive constant, R c and A p are constant, and q > 1 
(g 7 ^ 2). Inflationary models in such a power-law type gravity with q < 2 has also been 


examined in Ref. 


651 ]. Through the same procedures as those executed for the previous case 


in Eqs. (111.41) and (III. 5[) , namely, by deriving the Hubble parameter at the inflationary stage, 
the scale factor, the solution of the inflaton ip, the slow-roll parameters, and the number of 
e-folds, if n is close to 2 , we find 


A 


c- 


0?-i) s 


n R 


167 r 2 g (2 — q) 

8(2-9) 


pi exp 


2 -Np-^‘ 


(9-i r 


h) 




16(2 -qf 


(11.9) 


3 (g — 1) ’ 3 (g — l ) 2 

Indeed, for q = 1.99, we find n s = 0.962 and r = 1.08 x 10” 3 . Therefore, this inflationary 
model can be compatible with the Planck analysis. 


D. Reconstruction method of F(R ) gravity in the Jordan frame 

In this subsection, we reconstruct the form of F(R) in the Jordan frame. We note that 
cosmology in the Einstein frame may differ from that in the Jordan frame due to their 
physical non-equivalence. Hence, it is more convenient to consider these theories in the 
Einstein and Jordan frames as different cosmological theories. Here, we discuss inflation 
in F(R) gravity without its transformation to scalar-tensor theory. Eventually, the results 
may be different. Nevertheless, we demonstrate that F(R) inflation is also consistent with 
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the observations by the Planck satellite. From the other point of view, for the fairness, it 
should also be remarked that there are the debates on the issue of the equivalence between 
the (Jordan and Einstein) conformal frames in Refs. 66^ 68]. Especially, the investigations 
68 | seems to support equivalence of conformal frames for inflationary scenarios. It 


in Ref. 


is an issue of presentation rather than substance, but in the interest of fairness other points 
of view should be mentioned. 


The reconstruction method of F(R) gravity pro pos ed in Ref. 69J is as follows (for another 
reconstruction method of F(R) gravity, see Refs. 70lJ73|). We consider the action of F(R) 
gravity with matter action Setter as S — f d A Xy/—g [F(R )/ (2k 2 )] + Scatter- We define the 
number of e-folds as N = In (a*/a) with a* the scale factor at the fiducial time t*. We 
define G(N) = H 2 (N), so that R can be expressed as R(N) = 3 [(dG(N)/dN} + 4(5(IV)]. 
By solving this equation inversely, we get N = N(R). In the flat FLRW space-time, the 
Friedmann equation can represented as the second order differential equation of F(R) with 
respect to R, given by 


-9 G(N(R)) (4Gx(N(R)) + G m (N(R))) F rr (R ) 
f 3 (g(N(R)) + ±Gfi{N(R)) y ) Fr(R) - P matter = 0 , 


(IL10) 


where F RR = d 2 F(R)/dR 2 , G N = dG(N(R))/dN, G m = d 2 G(N(R))/dN 2 , and p matter is 
the energy density of matter. 


As an example found in Ref. 


54j . we study an exponential form G$ = H 2 (N ) = G 1 e' rAr + 


G 2 , where G\ (< 0), G 2 (> 0), and r (> 0) are constants. For this expression, we have 
e rA — [R— YIG 2 ) / [3Gi (4 + t)] . When the matter contribution is negligible, namely, 
Pmatter = 0, the solution of Eq. (III. 10]) is derived as 


F(R) = QiF (u + , w_, Z; t9) + Q 2 (12 G 2 - R) 


(1+1/r) F (1 + o;_ + -, 1 + o;+ + -, 2 + id 

T T T 


Here, uj± and § are defined as 

-3t — 2 ± Vr 2 - 20r + 4 


u± = 


4r 


d = 


12 G 2 - R 
3 G 2 (4 + r) ’ 


( 11 . 11 ) 


( 11 . 12 ) 


with F(qi,^ 2 , ^ 3 ] id) the hypergeometric function, where q (i = 1,...,3) are constants. If 
(N, Gi, G 2 ) = (50.0,-1.10,10.0) and (60.0,-1.20,15.0), we obtain (n s ,r) = (0.963,6.89 x 
















10 2 ) and (0.965,5.84 x 10 2 ), respectively 2 . Therefore, this model can yield the values of 
n s and r indicated by the Planck analysis. 


III. TRACE-ANOMALY DRIVEN INFLATION IN MODIFIED GRAVITY 


In this section, we review inflation by the quantum anomaly in the framework of F(R) 


gravity by following Ref. 


27]. The effect of the trace anomaly on inflation in F(T) gravity 




with T the torsion scalar in teleparallelism has also been studied in Ref. m (the explanati 


ions 


of teleparallelism exist, e.g., in Refs. [15 


, 120 )- 


A. Quantum anomaly 


It is known that the quantum anomaly appears via the procedure of the renormalization. 


For four-dimensional space-time, the trace o 


from the quantum anomaly becomes 


75 


79] 


the energy momentum tensor T^ A ' ) originating 


< T ^ a)11 ) = «i (Vv + 


a 2 G + a 3 DR , 


y\; — npvpan _ JDpvpcr p 
v y — ^ fits per 11 1 L pvpcr 

g = R^R pvpa - 4 RTR^ + R 2 


2 R7 v R pv + -R 


(HI-1) 

(HI-2) 

(III.3) 


Here, the brackets () denotes the vacuum expectation value. Moreover, is the Rie- 

mann tensor, R^ v is the Ricci tensor, R is the scalar curvature, C^ pa is the Weyl tensor, to 
whose square W corresponds, Q is the Gauss-Bonnet invariant, and □ = .g / " y V // V 1 , with V M 
the covariant derivative associated with the metric tensor g pv is the covariant d’Alembertian. 
In addition, the coefficients are defined as an = (Ns + 6 Np + !2Ny) / (19207T 2 ), a 2 = 
(Ns + HYf + 62iY v ) / (57607T 2 ), a 3 = —Ay/ (967T 2 ) with the number of real scalar fields 
Ns, that of the Dirac (fermion) fields Np, and that of vector fields Ay, where we have 
neglected the contributions from gravitons and higher-derivative conformal scalars. 

We set the values of aq and a 2 positive, but the qualitative consequences do not depend 
on aq, a 2 , and aq. For example, in the Af = 4 SU(N) super Yang-Mills theory, we have aq = 

2 The running of the spectral index a s = dn s /d In k is also estimated as a s = —5.06 x 10 -5 and —4.51 x 10 -5 
for (TV, Gi,G 2 ) = (50.0,-1.10,10.0) and (60.0,—1.20,15.0), respectively. 
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«2 = N 2 / (647r 2 ) (> 0) and 03 = —N 2 / ( 967 r 2 ), where we have used Ns = 6 N 2 , Np = 2N 2 , 
and Ny = N 2 with N 3 > 1. Here, (2/3) aq + a^ = 0. However, if the action has an additional 


R 2 term as 


80] [oqfV 2 / (1927r 2 )] J d 4 Xy/—g R 2 , where a 4 (> 0) constant, we 

hnd (2/3) aq + a 2 = — aqfV 2 / ( 167 T 2 ). In the classical level, the vacuum expectation value 
of ^ T$ A) ^ in Eq. (IIII. ID can be regarded as a contribution of matter in the right-hand 
side as R, IIU — (1/2) g IJiU R, = k 2 ^ Its trace reads R = a 1 [W + (2/3) HR] — a 2 Q + 

a^lUR + \a^N 2 K 2 / ( 167 T 2 )] QR. Accordingly, for the Yang-Mills theory in the curved space- 
time, the R 2 term plays a role of correction of the higher curvature to the Einstein gravity 
or it contributes to the energy-momentum tensor as matter. 


B. F(R) gravity with the quantum anomaly 


The action describing F(R) gravity is given by 


s = /<*V=9^r + s (Qi >, 


2 ( ot&N 2 \ 2 


(III.4) 

(HI-5) 


F(R) = R + 2k j I^^\ R‘ + /(«). 

where S (QA > is the action of the quantum anomaly and f(R ) is a function of R. For the 
original Starobinsky inflation 21], f(R) = 0. From the action in Eq. ( IIII. ip . the gravitational 
field equation reads 


R 


-RRfi v + —R 2 g fJ , l/ + V U R — g^OR 2 


fn(R) (r^ - \Rg»)j + \g»v ( f(R ) - Rfn{R)) 

(V M V, - g^U) f R (R ), 


(III. 6 ) 


where f R (R) = df(R)/dR. The trace of this equation is expressed as R = 
-k 2 {cqW - a 2 g - [cmN 2 / (16tt 2 )] DR} - 2 f{R) + Rf R (R) + 3Df R {R). 

In the FLRW background, the gravitational field equations are represented as 


4 H 2 = Pe S , ~^(2H + 3H 2 ) =P C 


Is- (HI-7) 

K* K* \ / 

Here, p e ff and P e ff are the effective energy density and pressure and they obey the equation 
of the conservation low p e g- + 3 H (p e ff + P e ff) = 0. Their expressions are given by 


P,<s = P (QA) + ^5 ( Rfa(R) - f(R) - 6 H 2 { r (R) - 6 Hfn(R) ) , 


(III.S) 
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P '» = -P (QA> + ^2 


-FtfiAFt! + f(R) + (4 H + 6 H 2 )f R (R) + 4 Hf R (R) + 2 f R (R) , (III.9) 


with the contributions from the quantum anomaly to the effective energy density and pres¬ 


sure 


'< Q M+ fa 2 p4 -®)( 


p(QA) _ _P_ 
- 3a 4 


a 2 ( 6 H 4 + 8 H 2 H 


18 H 2 H + 6 HH - 3H- 

2 ' 


f a A N‘ 
l 167T 2 


9 H 2 + 12 HH + 2 H+ 18 H 2 H 


(111.10) 

( 111 . 11 ) 


where p is an integration constant, and in deriving these expressions, we have used the 
conservation equation for p e g and P e fr- The p term corresponds to the energy density of ra¬ 
diation of the quantum state 80]. In what follows, we take p = 0 because at the inflationary 
stage around the Planck scale, the contribution of radiation can be neglected in comparison 
with that of the quantum anomaly as well as that of deviation of modified gravity from 
general relativity. 


C. de Sitter solutions by the trace anomaly 


We consider the case that f(R ) is given by an exponential form [49|, [50|, [81] 


f(R) = -2A C 


i / R 

1 — exp ( - — 

Ihc 


(III.12) 


where A c (> 0) and R c (> 0) are positive constants. For R/R c <C 1, i.e., in the late-time 
(e.g., present) universe, f{R) approaches zero, and therefore our model becomes equivalent 
to R 2 gravity with the quantum anomaly. While for R/R c ^ 1, namely, in the early universe 
such as the inflationary stage, the term of A c plays a role of the cosmological constant. When 
we expand the exponential term as exp (—R/R c ) = 1 — R/R c + 0((R/R c ) 2 ), these terms 
make the de Sitter solution realized by the trace anomaly unstable. This property can lead 
to the graceful exit from inflation. 

We derive the de Sitter solution at the inflationary stage. In the limit R/R c ^ 1, from 
Eq. (IIII.12I) we get f(R) ~ — 2A C . For this limit, the Friedmann (first) equation in (1III. 71) 
becomes (3 /k 2 ) H 2 k, 6 a 2 H 4 — 3a^N 2 / (167T 2 ) ^6 H 2 H + 2 HH — H 2 j + (A c /k 2 ). Solving this 
equation, we acquire the de Sitter solution 


H, 


de Sitter 


l 

4«2 k2 



8A c q;2^ 2 

3 


(III.13) 
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where we impose the condition A c < 3/ (8a 2 n 2 ) with «2 > 0 so that the solution can be real. 

We examine the instability condition of the de Sitter solution. If the de Sit¬ 
ter solution describes inflation, it should be unstable because inflation has to end. 
We represent the perturbation as H = H deS[ttei + SH(t), where \SH(t)\ <C 1. By 
combining it with the Friedmann equation, we obtain SH(t) + 3H de g it ter SH(t) = 
— [a 4 N 2 / (167T 2 )] [(1/k 2 ) — 4a 2 H^ e Sitter ] 5H(t), where we have neglected the terms propor¬ 
tional to exp (—R/R c ) in Eq. (IIII.121) because the stability of the solution is only related to 
A c . The solution for 5H(t) is written as 5H(t ) = Hex p (A ±t), where H is a constant and 
A± = 3H de s itter ± y/D^j /2 (the subscriptions ± of A± correspond to the signs of ± in the 
r.li.s.). The de Sitter solutions are unstable (and adopted to describe the inflation) only if 
the value of A + is a real and positive number, namely, 

9 64vr 2 1 9 

D — 9 H de Si tt er — ^ jy 2 J > 0 , J = — — Aa 2 H de Sitter > 0 ) (III. 14) 

where a 2 > 0 and a 4 > 0 have been used. 


D. Trace-anomaly driven inflation 


We investigate the observable quantities of the spectral index n s of the power spectrum 
for the the scalar mode of the density perturbations and the tensor-to-scalar ratio r in the 
trace-anomaly driven inflation in exponential gravity, namely, inflation is described by the 
de Sitter solutions in Eq. (IIII.13D . In the slow-roll inflation, for the exponential form of f(R) 
in flIII.12jl . we have 


tr e “A c a 2 K 2 (u + 2) 


<1, \r)\ 


u 


2 Ah 


< 1 , 


(III.15) 


AT2 (1 - 8A c a 2 « 2 /3) 

where 1 <C u < (l + y/1 — y) / (l — y/1 — y) with y = 8A c a 2 ^ 2 /3. As a result, we acquire 

1 


A 2 = 

n 32vr 2 a 2 e 


1 + \ 1 - 


8A c a 2 ^ 2 


u 


6 u 2 e u A c a 2 n 2 [u + 2) 


ns 1 N e N 2 (1 - 8A c a 2 « 2 /3) ’ 


r = 


16 u 2 e u A c a 2 K 2 (u + 2) 


(III.16) 

(III. 17) 
(III.18) 


N* (1 - 8A c a 2 « 2 /3) ' 

For u — 3, A c a 2 K 2 = 0.125, and N e = 76, we acquire n s = 0.960 and r = 1.20 x 10~ 3 . 
Consequently, the trace-anomaly driven inflation in exponential gravity can explain the 
Planck results. 
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IV. R 2 GRAVITY IN LOOP QUANTUM COSMOLOGY (LQC) 


In this section, we review R 2 gravity and its cosmological dynamics in LQC with the 


holonomy corrections along the investigations in Ref. 


A. F(R) gravity in LQC 


We explain F(R) gravity in the framework of LQC 87H89j. We consider the Einstein 


frame as in usual LQC only for the FLRW background with its spatially flatness [90]. In 
this case, the relation {/3lqc, Woiume} = 7bi/2 is satisfied 9l|. Here, { } denotes the Poisson 
bracket in terms of classical variables /3lqc = 7bi H with yei the Barbero-Immirzi parameter 
and the volume Koiume = a 3 , where a = x/Ero. Moreover, /3 lqc and Woiume are canonically 
conjugated quantities with each other, and these variables are the unique combination for a 
loop quantization |q 2|. We note that the hat shows the quantities in the Einstein frame. 

It is necessary to take the Hilbert space, in which the quantum states are described 
by (almost) periodic functions, so that the property of the discrete space can be includec 


35 


93 


For this purpose, we use the Hamiltonian with the general holonomy corrections 
Namely, we introduce A LQC = \J 7bi and replace /3 LQC with sin (a L q C /3 L qc) /A L qc 
in the Hamiltonian H = -3 (pIqc/'Ybi) Koiume + Kphime [(1/2) (dfi/dt) 2 + V($)\ = 0 with 
dt = y/F^dt by retaining {/3lqc> Woiume} = 7bi/2 [94j-l96|. As a result, we acquire the novel 


Hamiltonian "Hlqc- With the Hamiltonian equation dV vo \ ume /dt = {Woiume,'H lqc} and the 
Hamiltonian constraint 'Hlqc = 0, the Friedmann equation with the holonomy corrections 
reads [91 ] 

u 2 = \p i 1 - R—) ■ (tvi) 

<3 \ Pcritical / 

Here, p is the energy density of matter, and p C riticai = 3/ (Alqc7bi) 2 is the critical energy 
density. 


3 For LQC in teleparallelism, finite-time future singularities 
bouncing behaviors have also been studied in Refs. 


85 


82 


have been examined in Ref. [8j|, and 
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B. R 2 gravity in LQC 


For R 2 gravity, whose action is given by S — j d 4 Xy/—g [F(R)/ (2k 2 )] with F(R) = 
R + as r 2 R 2 , there appears curvature singularities in the early universe. In what follows, 
when we consider R 2 gravity, we analyze this action. On the other hand, for R 2 gravity in 
the context of LQC, it is possible that no singularity happens. We show this point below. 
It follows from Eq. (IIV.1D that 0 < p < p critical and - VPcriticai/12 < H < ^critical/12. 
In addition, for R 2 gravity, we have V(ip) = [1/ (8«sk 2 )] ^1 — exp ^—^/ 2 / 3 Kt/j j in the 
Einstein frame. Therefore, we find 0 < ( ddp/diy < 2p critica i and 0 < V(p) < /critical- For 
LQC, the Raychaudhuri equation becomes dH/dt = — (1/2) ( d(p/dt ) 2 (1 — 2/3/ /O cri ticai) - With 
this equation, we hnd dH/dt < (1/2) ( dip/dt)~ < /critical, and thus, it follows from this 


relation that 


R 


— ^/critical- 

Furthermore, the potential V((p) in the Einstein frame described above obeys 
dV((p)/dip = (1/Fr) \/V((p)/ (3ctgK 2 ). Substituting this equation into the relation be¬ 


tween R and R as R = F, 


R 


R + ( ddp/diy + a/6 (dV((p)/dLp) 


and using the relation 


1 — 2ctgK 2 


R + ( dtp/dt ) 

1*1 < 


> 1 - 18a s « 2 Pcriticai, we acquire 


18/critical d 


18 a g K 2 /critical 


2/critical \ 


asR 


- 


(IV.2) 


For a gK 2 <1/ (18/ cri ticai), the absolute value of R is bounded as in Eq. (1 1\ r .21) . In addition, 
the relation between H and H is represented as H = \/Fr H — (l/\/6) ( ddp/dt ) . Accord¬ 


ingly, \H\ is bounded. As a result, 


H 


= (1/6) | .R — 12H 2 1 does not diverge. Thus, there 


does not appear any singularity in R 2 gravity for LQC. 


C. Loop quantum R 2 gravity in the Einstein frame 


We further analyze the dynamics of R 2 gravity (i.e, F(R) = R + as r 2 R 2 ) in the 

Einstein frame, where equations become simpler than those in the Jordan frame. The 

equation of motion for (p is expressed by d 2 (p/dt 2 + 3 Hyddp/dt) + dV(<p)/dip = 0 with 

V (</) = [1/ (8«gK 2 )] ^1 — exp ^2/3. We introduce a variable T = exp (^/2/3. 

From the equation of motion for </, we obtain 

2 


d 2 T 

dt 2 


-T 



+ 3H —-T 
dt 


6ctgK 2 


T — 1 = 0 . 


(IV.3) 
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Since this equation is invariant under the transformation (H,t) (— H,—t ), the so¬ 

lution orbit draws the symmetric trajectories of the expansion and contraction phases 
on the ('&, dty/dt) plane in terms of the dA>/dt = 0 axis. That is, if there is a 
trajectory ('F(f), d$>{t)/dt) for the contraction phase (H < 0), we have a trajectory 
(T(—£), — rfT(— t)/dt) for the expansion phase (H > 0). 

In addition, the energy density is expressed as 

" / A \ 2 


P 


3 

4T 2 



+ 


1 

6agK 2 



This implies that H = 0 at (T(f), d^>(i)/dt) = (1,0). For p = p C riticai, we find 


4pcritical/ (3A) 


+ 


(■i-i/h ) 2 

80 s^^Pcritical/A 2 


1, 


(IV.4) 


(IV.5) 


with A = 1 — 8asK 2 /} critica i. This depicts an ellipse for A > 0, a parabola for A — 0, and 
a hyperbola for A < 0. There exists only the critical point at (T(f), dA>(t)/dt) = (1,0), 
where p = 0. All of the trajectories start from this point and come back to it. Therefore, it 
corresponds to both the beginning and end points of the universe. 

As a consequence, thanks to the holonomy corrections, in the Einstein frame, the bounces 
can occur when p = ^critical- The universe evolves from the contraction phase (H < 0). Its 
trajectory oscillates near the critical point (\l/(f), d\l/(t)/df) = (1,0) and the oscillatory 
amplitude becomes large. Eventually, the trajectory approaches the line p = /^critical and the 
bounce happens. After that, the expansion phase (H > 0) begins and the trajectory goes 
back to the critical point (T(f), dT(t)/df) = (1,0) with its oscillating behavior. 

Next, we explore the slow-roll inflation. With the slow-roll approximations ( dp/dty <C 
V(p) and \d 2 p/dt 2 \ <C 3 H [d(p/dt) , the Friedmann equation with the holonomy correc¬ 
tions and equation of motion for (p are written as H 2 = (1/3) V(p) (1 — E((p)//5 cri ti C ai) and 
3 H ( dp/dt ) + dV(<p)/dp = 0, respectively. The number of e-folds during inflation is derived 
by N e = 4> Hdi = f* H/ (dp/dt) dp « Jj‘ [V(0)l (dV(p)/dp)] (1 - V(0)/ftrttai) dp. For 
V(p) = [1/ (8as« 2 )] fl — exp \ , we find N e . ~ (3/4)exp (\/2/3n(p\] . More- 


*3Kipjj , we find N e w (3/4) exp 
over, for iV e >l, the slow-roll parameters are described as 


3 [1-1/ (4a S K 2 p cr iti C ai)] , __1_1_ 

41V 2 [1 - 1/ (8a S « 2 p crit ical)] 2 ’ V ~ iVe [1 — 1/ (8«S« 2 Pcritical)] ' 


(IV.6) 
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By using these expressions, the spectral index of scalar mode of the density perturbations 
n s — 1 = —6rj + 2 fj and the tensor-to-scalar ratio f = 16e become 


N e [1 - 1/ (8a S K 2 Pcritical)] ’ 

12 [1-1/ (4q S ^critical)] 

N 2 [1 - 1/ (8a S K 2 Pcritical)] 2 ' 


(IV.7) 

(IV.8) 


In the limit p C riticai —> oo, these expressions of h s and f become equivalent to those for pure R 2 
gravity without the holonomy corrections, i.e., the original Starobinsky inflation model. For 
instance, if N e = 68.0 and 8as^ 2 Pcriticai = 8.50, we acquire n s = 0.967 and f = 2.55 x 10“ 3 . 
Thus, in an R 2 gravity model with the holonomy corrections in the context of LQC, the 
spectral index of scalar mode of the density perturbations and the tensor-to-scalar ratio can 
be compatible with the Planck data. We mention that in the Starobinsky inflation model, 
for N e = 60.0 (68.0), we have h s = 0.967 (0.971) and r = 3.33 x 10" 3 (2.60 x 10" 3 ). 


D. Loop quantum R 2 gravity in the Jordan frame 


Furthermore, in the Jordan frame, we analyze the cosmological behaviors in R 2 grav¬ 
ity for LQC. For the comparison with the consequences in the Einstein frame, the point 
on the (T, d^/dt) plane at which the bounce occurs, i.e., H becomes zero. From H = 

Vf^\h 


(l/y/6) ( dRp/dt ) , if H = 0, we get H 2 = (1/4) fd&/dt) /T 2 . This equation 


leads to 


+ 


T - C±/B± 


= 1 


- -, (IV.9) 

/^critical/ (12U-I-) 2Ctg K~ /^critical / 

with B± = 1 ± A/8asft 2 p cr iticai anc ^ C± = 1 ± A/2asK 2 p cr i t i C ai, where the subscription ± in B± 
and C± corresponds to the sign ±. The case of + sign is for 0 < T < 1. In this case, this 
curve draws an ellipse for B + > o, a parabola for B + = 0, and an hyperbola for B + < 0. 
On the other hand, the case of — sign is for T > 1. In this case, it shows an ellipse. If 
the trajectory intersects this curve in the Einstein frame, the bounce happens in the Jordan 
frame. At the bouncing point, the relation H = (1/2) (d^/dtj /T has to be met. 

In the Einstein frame, the universe first contracts and finally expands at the critical point 
(T(t), d^>(t)/dt) = (1,0). With the relation between H and H and its time derivative 


H = ^77 



(IV.10) 


16 











(IV.ll) 


H = 


1 dA> ( ~ 1 d4//dt\ 


H 


J 


+ T 


dH 1 d fdA>/dt\ 

dt 2 dt\ T J 


2 dt \ 2 vp 

it is seen that in the Jordan frame, the universe begins and ends at the point (77, H) = (0, 0). 
It should be emphasized that the holonomy corrections yield the bounce in the Jordan frame, 
and hence, if they are absent, a singularity appears at the early stage of the universe. 


V. BOUNCING COSMOLOGY IN F(R) GRAVITY 


In this section, we review the cosmological bounce from F(R) gravity. Especially, we 


present the consequences found in Refs. 


41 


42], The bo uncin g behaviors in various modified 


gravity theories have also been investigated in Refs. (971-102]. We show that it is possible to 


reconstruct an F(R) gravity theory in which the matter bounce can happen in the framework 
of LQC. 

For the Friedmann equation with the holonomy corrections (1IV.ll) . the energy density of 
matter can be represented as p — p m / [(3/4 t 2 + 1] with p m a constant. In this case, the 
scale factor and the Hubble parameter read 86, Q 


a(t) = ( -p m t z + 1 


1/3 


m = 


(1/2 )p m t 


(V.l) 


(3/4) p m t 2 + 1 ' 

For these expressions, the solution of Eq. (III. 101) yields the form of F(R) to realize the matter 
bounce described by a and H in (IV.ll) for LQC. We solve Eq. (IV.ll) . and consequently acquire 

F(R) =1 1 R + 1 2 R 1/2 , (V.2) 

where X\ and X 2 are constants. We can set X\ = 1, so that the Einstein-Hilbert term should 
be included. 

In Ref. 42], the reconstruction of various modified gravity theories including F(R), F(Q), 


and F(T) gravity theories has been performed, where F(Q) is an arbitrary function of the 
Gauss- Bonn et invariant Q, to describe the two-times bouncing phenomena, called super- 
bounce 


103 


104]. and the ekpyrotic scenario |l05j in the context of LQC. 


VI. CONCLUSIONS 


In the present paper, we have reviewed inflationary models in modified gravity theories 
such as F(R) gravity including R 2 gravity with extended terms so that we can generalize 
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the Starobinsky inflation in R 2 gravity and derive its important properties to be useful clues 
to obtain the information on physics in the early universe. 

First, we have studied inflationary cosmology by modification terms of gravity, especially, 
inflation in F(R) gravity, by following Ref. 26]. The Starobinsky inflation in R 2 gravity is 
considered to be the seminal and significant idea of inflationary models in modified gravity 
theories. We have made the coformal transformation from the Jordan frame (namely, F(R) 
gravity) to the Einstein frame (i.e., general gravity plus the scalar field theory), and given 
slow-roll dynamics of inflation in the Einstein frame. In addition, we have reconstructed 
F(R) gravity models, which are an extended version of the Starobinsky inflation model in 
R 2 gravity and general relativity with power-law correction terms. 

Second, we have explored the trace-anomaly driven inflation in F(R) gravity along the 


discussions in Ref. 


27]. We have first explained the quantum anomaly appearing through 


the process of the renormalization in four-dimensional space-time. We have further discussed 
F(R) gravity with the quantum anomaly and the de Sitter solutions for inflation due to the 
trace anomaly. 

Third, based on Ref. { 28 ], we have examined inflation in R 2 gravity and the cosmological 
evolutions for LQC with the holonomy corrections. We have analyzed R 2 gravity for LQC 
in both the Einstein and Jordan frames. We have found that in the Jordan frame, owing to 
the holonomy corrections, the bounce can happen, and accordingly the cosmic singularities 
can be removed, although such singularities appear in ordinary R 2 gravity. 

In these three inflationary models, we have shown that the spectral index of scalar modes 
of the density perturbations and the tensor-to-scalar ratio can be compatible with the Planck 
analysis. 

Furthermore, we have presented the recent developments of the bounce cosmology in 
F(R) gravity obtained in Refs. 41, Q. It has been performed that an F(R) gravity theory 
can be reconstructed, where the matter bounce occurs in the context of LQC, thanks to the 
reconstruction method of F(R) gravity. Recently, the reconstruction of F(R), F(Q), and 
F(T) gravity theories have also been executed, in which the super-bounce (i.e., two-times 
bounce behaviors) and the ekpyrotic scenario for LQC can be realized. 

In this work, we have concentrated on the accelerating universe from F(R) gravity. Note, 
however, that it is possible to extend this study for more complicated versions of effective 
gravity, which comes from quantum gravity. Particularly, it has recently been demonstrated 
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that successful inflation consistent with the Planck data may e merge from multiplicatively 


renormalizable higher derivative quantum gravity in Ref. 


106], 
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